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ABSTRACT 


In this thesis, a Linear Quadratic Gaussian Controller (LQG) is designed for a tail 
controlled surface-to-air missile model in order to meet design specifications. The 
mathematical model of the flexible missile is subject to uncertainties that may arise from 
unmodelled dynamics, parameter variation or linearization of nonlinear elements. Since 
these uncertainties are not taken into account in the LQG controller, ,: Analysis is 
applied to the design in order to evaluate the Robust Performance, Robust Stability, and 
Nominal Performance of the system. Finally, a Linear Quadratic Gaussian controller is 
designed using Loop Transfer Recovery (LQGLTR) in order to improve the Robust 
Stability of the system. It is found that the Robust Stability of the design is improved, but 
as a consequence of losing nominal performance. The p Analysis and Synthesis Toolbox 
and the Control Toolbox of MATLAB were used for the design, assembly, analysis and 


simulation of the missile flight control system. 
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I. INTRODUCTION 


Highly evasive target maneuver levels and increasingly effective electronic 
countermeasures push current air defense missile systems to the limits, resulting in more 
stringent autopilot performance requirements. Conventional autopilot design techniques 
have worked well in the past, but new design methods are required to obtain improved 
performance and robustness characteristics from the flight control system in order to 
satisfy future design specifications. A significant research effort has been directed toward 
the design and implementation of robust controllers which can guarantee both stability 
and performance 

Robust control design methods such as H, Optimal control or Linear Quadratic 
Gaussian Loop Transfer Recovery optimize performance and stability based on 
engineering models which include performance specifications and descriptions of how 
uncertainty modifies the nominal plant. 

Recent advances in robust feedback control, in particular H, contro! and 
-synthesis, allow for new design methodologies to be applied to the design of missile 
flight control systems. H, optimal control provides the basis for controller synthesis 
while mu analysis characterizes performance and stability in the presence of a defined 
structure for uncertainties. Recent advances in robust control theory offer good prospects 


for meeting the design needs of next generation missiles. Several anticipated benefits of 








the robust control design appioach are: greater flexibility in the choice of airframe 
geometry, full use of available airframe maneuver capability and greater tolerance to 
uncertainty ia design models. 

Linear quadratic controllers using state estimate feedback are optimal for the 
nominal model but the performance may be far from satisactory in a real life situation in 
which the plant differs somewhat from the model. Attractive passband robustness 
properties of full-state feedback optimal quadratic designs may disappear with the 
introduction of a_ state estimator. The  Linear-Quadratic-Gaussian with 
Loop-Transfer-Recovery (LQG/LTR) methodology provides an integrated frequency 
domain and state approach for design of multi-input multi-output (MIMO) control 
systems. The advantages of the methodology lie in its ability to directly address design 
issues such stability robustness and the trade-off between performance and allowable 
control power. 

In this study the problem of the design of a pitch plane autopilot to track normal 
accelerations commanded from the guidance system for a tail controlled flexible missile 
[Ref. 1] using Linear Quadratic Gaussian techniques is presented. 

In general when a missile is flying at an angle of attack alpha, lift is developed. 
This lift may be represented as acting at the center of pressure of the structure. The missile 
will be statically stable or unstable without corrective tail deflections depending on the 


location of the center of pressure relative to the center of mass. 


The control problem requires that the autopilot generate the required tail-deflection 


delta to produce an angle of attack, corresponding to a maneuver called for by the 
guidance law, while stabilizing the airframe rotational motion. The primary objective of 
the missile autopilot is to control the missile's flight path while meeting the specified 
performance requirements and to mantain the specified stability margins of the airframe 
over the entire flight envelope, all in the presence of parameter variations, unmodeled 
dynamics and disturbances. The missile flight control system must guarantee stability 
and performance in the face of large aerodynamic uncertainties and disturbances. This 
requires the feedback controller to maintain system stability and loop performance for all 
variations in the plant behavior. The autopilot design goals are simply to make the missile 
respond as fast as possible while maintaining stability. Stability robustness goals include 
robustness to neglected bending dynamics and uncertainties in the aerodynamic 
parameters. 

A pitch plane model describing the longitudinal dynamics of the missile is available 
for this design. Reasonable accurate mathematical models of the rigid body transfer 
functions from tail-deflection to the sensor outputs are available too. Sensor 
measurements for feedback include missile rotational rates (from a rate gyro) and normal 
acceleration (from an accelerometer). 

The sensors are assumed to be subjected to noises, dynamic characteristics, 


calibration errors and drifts. 











The mathematical model of the flexible missile is subject to uncertainties thatmay 
arise from unmodeled dynamics, parameter variation, linearization of nonlinear elements. 
Model perturbations or uncertainties such as aerodynamic stability coefficient variation 
uncertainties and actuator perturbations are considered in the model. Multiplicative 
parameter uncertainties are considered upon the pitching moment and normal force 
aerodynamic stability derivatives with respect to angle of attack. 

In general there are other uncertainties associated with a missile system that can be 
considered but are not part of this specific model. For example, the missile may be 
constantly subjected to environmental disturbances such as wind gust. The wind may 
induce a sideslip excursion and affect the measurements. 

The problem of structural flexibility is also taken into account since this is a long 
and slender missile. This introduces problems relating to the elastic nature of the airframe. 
As the autopilot bandwidth increases, coupling between the flight control system and the 
flexible dynamics of the airframe may take place which limits the achievable airframe 
response and possibly destabilizes the system. Coupling occurs between the elastic modes 
and the control system as the control system gyros and accelerometers sense the flexure 
motion and the rigid body motion. The coupling arises from the fact that the attitude and 
rate gyros and accelerometers sense both the rigid body changes and the body bending 
motion. 

The controller must avoid saturating the tail deflection actuator rate capabilities and 


destabilizing high frequency flexible body modes of the missile. The push for a smaller 








time constant to achieve a rapid command response is in conflict with the requirement for 
robust stability against high frequency uncertain dynamics. During flight, aerolastic 
forces also act on the airframe, causing high-frequency vibratory modes. The vibrations 
affect both the acceleration and the pitch rate but only the second one is considered in the 
model of the missile. 

This work is organized as follows. Robust multivariable control theory as well as 
the Structured Singular Value will be introduced in Chapter II. The Structured Singular 
Value will play «- .mportant role in the evaluation of the missile's performance in 
Chapter IV. Linear Quadratic Gaussian techniques will be applied to the missile model in 
Chapter III to yield the desired specifications. This chapter also presents the missile and 
airframe models used in the design. A technique known as yp analysis is applied in Chapter 
IV in order to evaluate Robust Performance, Robust Stability and Nominal Performance 
of the Linear Quadratic Gaussian design for the model of the missile. Linear Quadratic 
Gaussian with Loop Transfer Recovery is applied in Chapter V to improve the stability 
robustness of the system. The final chapter summarizes the results and suggest further 


research. 














I. BACKGROUND 


A. GENERAL FEEDBACK PERTURBATION MODEL 
Any general feedback system in the presence of uncertainties can be represented by 
the general A-P-K framework shown in Figure 1, where the plant model P,,, has been 


expanded to include additional input and output variables. 


WwW P Y 
(s) 
K 
(s) 
P : Plant Model 
A : Uncertainty 
K : Feedback Compensator 
W_ : _ Disturbance Input 
Wd : Feedback through uncertainty 
Yd : Output Leading to feedback uncertainty 
Y : Reference output 


Figure 1. General A-P-K Framework. 




















The additional input Wd includes those plant variables which are subjected to plant 
parameter variations, while the additional output Yd denotes plant variables upon which 
those plant variations may act. The matrix A(s) represents a transfer function matrix 
describing the effect of plant uncertainties. Any linear time invariant contro! design 
problem can be arranged into this form [Ref. 2]. The convention adopted here is to 
normalize exogenous inputs W, outputs Y so that the uncertainty A is of a magnitude 
less than 1. This requires that all scalings be included in the plant description P. 

The plant is then modeled as having three inputs and three outputs and is 


completely described by the transfer function matrix. 





Yd Pir Pio Pig |} Wd 
Y {=| Par Pa Pas W i, (1) 
m P3; Pa. Pas 


or, in a compact form, 


Yd Wd 
Y {=P W I. (2) 
m U 


In a notation commonly used the generalized plant can be expressed as 
Py) =C(sl-A)"'D=| -.. + |, (3) 


where the matrices A, B, C, and D are given in the state space model of the plant [Ref. 1]. 








B. UNSTRUCTURED UNCERTAINTY MODEL 

A linear model of a multivariable linear system includes a description of the 
uncertainties and their structure. The uncertainties are due to parameter variations, 
unmodeled dynamics, neglected flexible body dynamics, non-linearities, and various 
assumptions that affect the model of the real system [Ref. 1]. 

In general, uncertainties in dynamic models could enter the system in a variety of 
ways, such as parametric or non-parametric, structured or unstructured, stochastic or 
deterministic. 

An unstructured uncertainty is modeled as an unknown with a bounded transfer 
function incorporated in series or parallel with the plant [Ref. 3]. This type of model of 
the plant uncertainty is termed unstructured since no detailed model of the perturbation is 
employed. Figure 2 shows the unstructured uncertainty models that will be considered. 

The actual plant will equal the nominal plant for each of these cases when the 
perturbation equals zero. For evaluation of stability, robustness or performance 


robustness these perturbations are assumed to be contained within bounds. 


||Ai(S)l] < Anax (4) 














Additive 


[e@}—S 


G(s) 


Input Multiplicative 





where 
G,(s) represents the nominal plant 
x(s) represents a perturbation 


Figure 2. Unstructured Uncertainty Models. 








C. STRUCTURED UNCERTAINTY MODEL 

Structured uncertainties refer to plants which contain one or more uncertain 
parameters, they also refer to plants which contain both uncertain parameters and 
unstructured uncertainties. 

Uncertain blocks and parameters from various plant locations can be rearranged as 


shown in Figure 3 [Ref,. 3]. 





Figure 3. Closed Loop System with Feedback Perturbation. 


In this figure M(s) is the combination of the nominal! plant P(s) with a feedback 
compensator K(s), and it represents the dynamics of the nominal closed loop system. It 


can be divided into 4 blocks as shown below [Ref. 1]. 


Ma Miz 
M= 5 
ae Maze () 


relative to the inputs and outputs 


10 


This matrix partition is called Linear Fractional Transformation and is derived from 


the following linear equations. 


(6) 


The term M,, may be viewed as the nominal closed loop dynamics, assuming no 


uncertainties and A as a linear-fractional uncertainty. The matricesM,,, M,,, M.,, F (M, 


14> 129 


A) describe how A affects the overall dynamics, expressed by the transfer function matrix 


F(M, A) = Moo + May AC — My, A) M ye (7) 





We assume independence of the terms in the structured perturbation which is of the 


form 


Ax(S) 0 0 
asy=| 9 428) -~ 0 (8) 
0 0 Ap(S) 
or, more compactly 
A=diag{A;, Az, ...Ap} (9) 


The individual uncertainties, A, that make up the A block may have one of several 


representations [Ref. 1]: 


1] 








1. Scalar Block; uncertainty appears only once in the real system, such as 
parameter uncertainty. 
2. Repeated Scalar Blocks; where the uncertainty appears multiple time in real 


system. 


3. Full Blocks; example of which includes multivariable neglected dynamics and 


the performance block. 
As in the standard model for unstructured uncertainties, the structured uncertainties 


will all be scaled so their infinity norms equal one [Ref. 3] 


WAll.. = rte =--- = [lApll. = 1. (10) 


where we define the infinity norm as||All,. =sup [AQjo)]. 


D. STRUCTURED SINGULAR VALUE 
1. Singular Value Decomposition 


Any m X n matrix can be factored as 
A=VirU' (11) 


where V is an m X m orthogonal matrix, U is ann x n orthogonal matrix, 2 is anm xn 


matrix of the special form 








where 


D= co (13) 


is diagonal. 
The elements 6, > 6, 2... 2.6, > 0 are all positive real, and are called the singular 
values of the matrix A. The singular decomposition has the property that the L; gain of 


the matrix lies between the smallest and the largest singular values. 


fIMx fl = 
G=Op < =] £01=6 





(14) 


where o is the smallest singular value and 6 is the largest singular value. 


For SISO systems the concept of "gain" is synonymous with the classical Bode 
magnitude plot. This Bode magnitude is equal to the amount of amplification of a 
sinusoidal input by the system. 

For multivariable systems, the extension of the Bode gain is the singular value of 
a transfer matrix M evaluated at a frequency upon the imaginary axis, M(jo) [Ref. 4]. 

If we consider the transfer function matrix M evaluated at any frequency point 
M(jo), the largest amplification that a sinusoidal signal of frequency will receive as it 


passes through M is equal to the maximum singular value of M at s = jo. 


13 





2. The Principal Gains 
The principal gains are the frequency dependent singular values of the system 
They can be used to assist the designer in determining the robustness of the multivariable 
system and are a good indicator of the system's ability to cope with disturbance rejection. 


The principle gains will be denoted as [Ref. 5] 
a(G(ja) = 7 (15) 


where G(jw) is the frequency dependent transfer function of the system. 
If the principal gains of the transfer function matrix of the MIMO system are 


ordered such that 


¥1 <Y2<...<Yn, (16) 


then y,and y, have the property that they bound theL, gain between the disturbance input 


vector d(s) and the output vector y(s) 


Thus if y, is determined for the transfer function matrix then this is the maximum 
gain which can be obtained by any disturbance vector d(s) and gives a worst case 
indication of the disturbance rejection properties of the control system while y, gives an 


indication of the best rejection by the system. 


14 














Large values of o,,, over a large range of frequencies result in good tracking of 
reference inputs and good rejection of disturbances. 


The maximum gain of the system will be denoted by 
o(G(jo)) (18) 
and the minimum gain by 
o(G(§ja)) | (19) 


The maximum gain that a stable linear system M can deliver is then the 
supremum of the maximum singular values taken over all frequencies. 
3. The Infinity Norm 


Mathematically, the infinity norm of a transfer function matrix M is defined as 
the supremum over all frequencies of the maximum singular valued of M(ja). 


It is denoted by 


IMI. =sup a[MGo)] (20) 


The first engineering relevance of the infinity norm is that if a control designer 
seeks to design a controller K for minimizing the transfer of energy or the transfer of 


power or the transfer of white noise through a closed loop systemM which depends on K, 
then the controller K should be designed so that the infinity normM, ||M|l,., is minimized 


over all controllers which stabilize the closed loop system M. 


15 











The second engineering relevance of the infinity norm is that it is essential in 
establishing robustness of uncertain systems. 
4. The Structured Singular Value 
The structured singular value of a complex matrix M denoted by p(M(ja)) , with 
respect to the block structured A, is based on the multivariable Nyquist Theorem and is 


defined as 


1 


H(MU0) = Trin fa): det + Mago) = 0 


(21) 


where A denotes the set of all transfer function matrices with the appropriate block 
diagram form. 

It can be shown that the structured singular value is the inverse of the stability 
margin for multivariable systems [Ref. 3]. In other words, M is the inverse of the smallest 
magnitude of a destabilizing perturbation of M. 

The structured singular value is a measure of robustness of the model, to 
structured perturbations. 

The structured singular value provides an indication of how much uncertainty 
can be tolerated before the systems becomes unstable. It will be used in Chapter IV in 
order to evaluate the robust performance and robust stability of the Linear Quadratic 


Guassian design for the missile model proposed in Chapter III. 


EE 


Ill. LINEAR QUADRATIC GAUSSIAN (LQG) DESIGN 


In this chapter we apply LQG techniques to the design of a controller for the model 
of a missile with flexible modes proposed in Section A, in order to meet design 


requirements specified in Section B. 


A. MISSILE MODEL 
The airframe model used in this thesis was presented by Bibel [Ref. 1]. A block 


diagram configuration for the model of the surface launched, antiair homing missile to be 





considered is shown in Figure 4. The missile is formed basically by the following 
sub-systems: fin actuator, rigid body dynamics, rate gyro, accelerometer, and flexible 
body dynamics. The State Space Model representation of this Multi-Input Multi-Output 


System can be described by 





x =Ax+Bu+IW, (22) 


y = Cx+Du+¥PV, (23) 





where 
A = Plant Matrix, B = Input Matrix, C = Output Matrix, D = Feedforward Matrix, 


IT = Input Matrix, ‘Y =Measurement Noise Input Matrix. 
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The reference input to the system 1s given by the commanded acceleration and the 
control input is given by the commanded fin deflection. Two disturbances entering the 
system, accelerometer bias and gyro drift, and twomeasurements, measured acceleration 
error and measured pitch, are modelled as white noise processes with intensities W and 
V, respectively. 


_| Vit) Viet 
vio=| V0) Volt | t2>to, (24) 


where V(t) and V,(t) are the intensity of W and V, respectively and are correlated. 
V,,(t) is equal to O (uncorrelated) and W and V are independent. 


We define the covariance matrices as 


E(W(Q)W/(t + 1)] = Q6(t- 1), (25) 


E[V(t)V/(t+2)] = R&(t- 1). (26) 


For the purpose of modeling these two disturbances, models for biases have been 
added to the nominal plant model. The error between the commanded acceleration and 


the measured normal acceleration and the measured body pitch rate are the feedback 


signals to the controller. 


The state variables considered in the nominal model are summarized in Table I. 
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TABLE I. NOMINAL MODEL STATE VARIABLES AND 
PARAMETERS. 


Fin Rate 
Actual Fin Deflection 
Angle of Attack 
Pitch Rate 
Commanded Fin Deflection 


Normal! Acceleration (measured by the 
accelerometer) 


Commanded Acceleration 


Error between n. and n 
Measured Acceleration 
Measured Pitch Rate 
Measured error in Acceleration 








The model to be considered in the design of the linear quadratic gaussian controller 
combines the nominal plant and the process-noise model but does not include the 
description of the uncertainties and their structure (i.e., two aerodynamic stability 
coefficient variation uncertainties, two actuator perturbations, multiplicative parameter 
uncertainties on the pitching moment, and normal force aerodynamic stability derivatives 


with respect to the angle of attack). 





B. AUTOPILOT DESIGN REQUIREMENTS 
The primary objective for the design of the autopilot is to obtain the fastest time 


response possible while maintaining stability robustness. The following time domain 
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specifications are to be considered in the design of the linear quadratic gaussian controller 


for the missile autopilot. 


Settling Time Ts < 0.4 secs (Time constant = 0.1 secs) 
Steady State Error e,,=0.1% (Due to a step command) 
Overshoot Ov < 4% (Due to a step command) 


The linear quadratic gaussian controller has to be able to handle the effects of the 


airframe's elastic dynamics. 


C. PROPOSED DESIGN 

The design of the linear quadratic gaussian controller for the model of this missile is 
given by a combination of the solutions of the Stochastic Optimal Regulator and 
Optimal Estimator problem. 

The optimal control, solution of the stochastic optimal regulator problem, which 


minimizes the next performance index 
J=JF y?(t) + pu2(tdt, (27) 


where y(t) represents the output of the accelerometer and u(t) the control input to the 


system is given by state feedback 


u = —K:x (28) 


K,=R“B/P (29) 
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where P is the unique positive definite matrix satisfying the Matrix Riccati Equation 


PA+A/P - PBR-'B/P+Q=0 (30) 


provided that Q is symmetric and positive semidefinite, R is symmetric and positive 


definite, the pair (A,B) is controllable. with the controllability matrix given by 


€ = [B AB A2B...A™"B] (31) 


The scalar p plays the role of a Lagrange multiplier. Asp decreases, the integrated 
square regulating error decreases, but the integrated square input increases. From 
simulation, p was chosen to be 0.1. The response of the missile due to a unit step 


commanded acceleration with the Linear Quadratic Regulator is shown in Figure 5. 


0 0.1 0.2 0.3 0.4 
Time (secs) 


Figure 5. Missile Response Due to a Unit Step Commanded Acceleration. 





The solution of the Optimum Observer or State Estimator given by Kalman and 


Bucy can be expressed by the following Differential Equation 




















x’ = Ax +Bu+K.(y - Cx). 


The Kalman Gain K, is optimally chosen so as to minimize 


Exe"(t)Wite()}. 


and is given by 


Ko(t) =-P. (CR (t): 


where P.(t) is the solution of the matrix Ricatti equation 


P.(t) = Pe(t)F/(t) + F()P.(t) - Pe(QH()R™ (QH/(t)P.(t) + Q(t), 


where 


e(t) = x(t) - X(t), 


and 


F.(t) = F(t) + K.(t)H/(t). 


This has a solution provided that 
W is symmetric and Positive Semidefinite 
V is symmetric and Positive Definite 
The pair (A,C) is observable 


where the observability matrix is given by 
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(32) 


(33) 


(34) 


(35) 


(36) 


(37) 








9=[CT ATCT A™CT..AT'CT] (38) 


The optimal output feedback control system is represented below. This combines 
the model of the nominal plant with the optimal feedback gain and the optimal estimator. 
It also includes a state command matrix N, that defines the desired value of the state, x. 
Since this is a tracking problem, N, should transform the reference input r to a reference 


state. It can be shown that N, can be calculated from the following matrix equation: 


iN LS ahi a (39) 


where ® and I represent the plant matrix and input matrix of the discrete representation 
of the model of the system. 


The Reference State is given by 
x,=N,¢ 


and it represents the state we want to track. 
Equations (22), (28), (32), and (40) can be combined in order to obtain the State 
Space Model for the Closed Loop Optimal Output Feedback System shown in Figure 6. 


From the Control Law given by 
u= —Ki(X —X,), (41) 


and substituting Equation (41) into Equation (22) we obtain the closed loop dynamics 








x = Ax - BK, (x - x,) + TW. 


ACCNOISE GYRONOISE 











A 
x =Ax+ Buel 
en a 


Missile Model 






R= Ax + Bu + K o(y-Cx) 






Optimal Kalman Filter 
Feedback 
Gain Be State Command 


Matrix 


Cmdacc 


Figure 6. Closed Loop Optimal Output Feedback System. 


which can be written as 


x = Ax - BK;x + BKx, + PW. 
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(42) 


(43) 











Combining these equations and performing the required algebra we obtain 
; r 
x! }-| ec A-(B nore -K.C [ 3 |.) (B oe 5 oe i (44) 
which represents the State Space Model of the Optimal Output Feedback Control! System. 
The acceleration time response of the closed-loop system to a unit step acceleration 
command was obtained by computer simulation and is shown in Figure 7. 


The missile's angle of attack, body pitch rate, gyro output, commanded fin 


deflection and fin rate, are shown in Figures 8, 9, 10, 11, and 12. 





“0 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 7. Achieved Normal Acceleration Due to a Unit Step Commanded Acceleration. 








0 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 8. Gyro Output Due to a Unit Step Commanded Acceleration. 





INT4 


0 0.2 0.4750 (sec) 0.6 





Figure 9. Pitch Rate Due to a Unit Step Commanded Acceleration. 











0.15 


0.1 


INT3 


0.05 


0 0.2 0.4_ 0.6 0.8 1 
Time (sec) 


Figure 10. Angle of Attack Due to a Unit Step Commanded Acceleration. 


INT2 
° 


0 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 11. Actual Fin Deflexion Due to a Unit Step Commanded Acceleration. 
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0 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 12. Fin Rate Due to a Unit Step Commanded Acceleration. 


The following time domain specifications were achieved with a standard LQG design: 


Final Value : 1.0000 unit 

Time toPeak : 0.4500 secs 

% Overshoot : 0.0000 % 

Rise Time : 0.2500 secs 

Settling Time : 0.4000 secs 
Steady State Error : 0.0000 units 


These specifications were obtained from Figure 7. The acceleration response is 
critically damped and is the fastest that can be obtained for this particular design. 

If we observe the missile's angle of attack and body pitch rate responses it can be 
noticed that these two do not affect the flexible body dynamics. The fin deflection limit 


of 40 degrees and the fin rate limit of 300 degrees/sec are not exceeded. 
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So far all of the required specifications for which the controller was designed have 
been met, but we still have to check the stability of the system. The LQG controller was 
designed for the nominal plant model of the missile and may or may not meet robust 
stability requirement. 

In the next chapter this design will be tested using a technique called p1-analysis 
which evaluates the robust performance, robust stability and nominal performance of the 


design. 
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IV. p-ANALYSIS FOR LINEAR QUADRATIC GAUSSIAN DESIGN 


The goal of any controller design is that the overall system is stable and satisfies 
desired performance requirements. These requirements should be satisfied at least when 
the controller is applied to the nominal plant, that is, nominal stability and nominal 
performance is required. In practice the real plant is not equal to the model. Even though 
there are uncertainties in the dynamics of any real system we try to quantify the 
unmodeled dynamics. 

In the design of the Linear Quadratic Gaussian controller, the nominal plant model 
of the flexible missile was assumed to be known exactly. This is rarely the case due to 
unmodeled dynamics, parameter variations, uncertainties, etc. Even if these uncertainties 
are neglegible, it may be desirable to design a controller for a plant whose parameters 
vary deterministically as a consequence of environmental factors or component failures. 

When the controller designed for a nominal plant model is implemented on the real 
system, there are no guarantees on the resulting performance of the system. Even 
requirements as basic as stability may not be met. The deviation from the expected 
behavior of the system clearly depends on the accuracy of the model. Because models are 
never perfect, robustness analysis is necessary to determine the possibility of instability or 
inadequate performance in the face of uncertainty of the plant dynamics. A controller 


should guarantee stability, proper response to commands, and reduction of response 
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perturbations caused by disturbance inputs. Ability to perform the first function under 
parameter variation is called stability robustness, while ability to conduct the remaining 


two functions is termed performance robustness. 


A. NOMINAL STABILITY AND ROBUST STABILITY 
If the LQG controller is included in the nominal plant P, the following class of 
general matrix transformation called Linear Fractional Transformation can be defined 


as [Ref. 4] 


FP, K) = P,, + Pi2K(l = Pook)! P., (45) 


Wd 
W 


F(P,K) denotes the transfer function rom| | to | . | This particular structure is 


shown in Figure 13. 


Yd Wd 





Figure 13. Linear Fractional Transformation. 
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The complex matrix M, represents the closed loop transfer function between 


far | and| Ms | It can be partitioned as 
Mis Mi2 
M= : 46 
Mor Moz “ 


Therefore, the system is nominally stable ifM is stable. The stability of a MIMO 
system is determined by analyzing the feedback system shown in Figure 14. The poles of 


this feedback system are the solutions of the charateristic equation: 


det{I + Mi (s)A(s)} =0 (47) 





Figure 14. Nominal Closed Loop Plant with a Feedback Perturbation. 


It can be shown, by the small gain theorem [Ref. 3], that the system is internally stable if 


and only if 
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[Mri HcWAll. <1 (48) 


for all possible perturbations. Given that 


All. $1, (49) 


the system satisfies the robust stability requirement if 


Marl. <1 (50) 
when A is a full uncertainty block or 
Maxp rc (M11 Gw))IAll.. < 1 (51) 
if 
AcA (Structured Uncertainty) (52) 


where jt is called the Structured Singular Value (SSV) defined in Chapter II. 


B. NOMINAL PERFORMANCE 
The closed-loop system achieves nominal performance if the performance objective 


is satisfied for the nominal plant model. Mathematically, the system satisfies nominal 


performance if 


Mell. < 1. (53) 








BR 


C. ROBUST PERFORMANCE 
A closed-loop system will be said to have robust performance if the total system 


meets some minimum objectives for all possible perturbations. Given a structured 








uncertainty 
; | 
AeA (54) | 
| 
where 
All. <1 (55) 
and the performance specification 
ITS). <1, (56) 
where 
y(s) = T(s)W(s), (57) 


the system will posses "Performance Robustness" if and only if 
Max p[M(jw)] < 1. (58) 


To compute the SSV, a block diagonal structured perturbation of the form| : 4 | 
Pp 


where AeA and A, is unstructured, is assumed. 











These are all criteria which quantify the performance and stability margins of a 
mulitvariable system. In particular, if we plot these as functions of frequency, the system 
meets all of the specifications, provided they are all less than one. The uncertainties are 
normalized as in Equation (49). All of this can be shown by using the missile model 
introduced in Chapter III with the LQG design. Robust performance, robust stability, and 
nominal performance tests are shown in Figure 15. Notice that only nominal performance 
is met, since it is less than one, but robust performance and robust stability are not met. 
This means that we can find a set of perturbations for which the system will not remain 


stable. 


Gain 





10° 10 10 10° 10° 10° 10° 


Freq (Hz) 


Figure 15. Robust Performance, Robust Stability, Nominal Performance of LQG Design. 
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The LQG controller meets the design goal at the nominal operating condition. 


However, it becomes sensitive in the presence of parameter variations. 


° 


Amplitude 


0.6 








“0 0.2 0.4 0.6 0.8 1 
Time (sec) 


Figure 16. LQG Design Under Uncertainty Variations. 





As stated in Chapter III, there are four model perturbations or uncertainties in the 
design model of the missile. These include two aerodynamic stability coefficient 
variation uncertainties and two actuator perturbations. One of the desired controller 
attributes is to remain insensitive to angle of attack changes during flight. Step responses 
of the closed-loop system with the LQG controller designed in Chapter III for different 


values of the uncertainty Z, are shown in Figure 16. It can be noticed that the LQG 
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design meets the specifications for which it was designed when the parameter Z, remains 
in its nominal value. Once the parameter Z, changes, the LQG controlller is unable to 
track to desired output, and the specifications are no longer met. 

The LQG regulator guarantees stability under nominal operating conditions and 
meets the design specifications for which it was designed. However, stability robustness 
was not considered explicitly in establishing the optimality criteria. When the estimator is 
part of the feedback loop of a stochastic regulator, dynamic interactions between the 
estimator and the controlled system increase the potential for instability. 

In the next chapter, a technique known as loop transfer recovery will be used in 


conjunction with the LQG design in order to improve the robust stability of the system. 
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V. LQG LOOP TRANSFER RECOVERY DESIGN (LQGLTR) 


A major objective of feedback system design is to achieve anominal performance 
specification for a given design model of the plant, and to maintain this performance over 
a range of expected errors between the design model and the true plant. 

As stated in Chapter IV, the attractive passband robustness properties of full-state 
feedback optimal quadratic designs might disappear with the introduction of a state 
estimator. However the loop properties of a LQ design can be recovered by a suitable 
adjustment to the LQG design process with a technique known as Loop Transfer 
recovery. 

Given the state space model of the system defined by the matricesA, B, C, and D, 


the transfer function of the system can be denoted by 


P(s) = COB +D, (59) 


where 


®=(sl-A)-'. (60) 


The open loop transfer function at the plant input can be written as 


L.(s) = KB (6 1) 
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where L,,, refers to the target loop, and K is the controller gain. By adding fictitious 
noise to the plant input model representing a loose way plant variations, uncertainty, or 
unmodeled dynamics, the estimator design can be adjusted in such a way that the 
robustness properties of the LQG design can be recovered. 

For this purpose the controller K(s) can be parameterized as a function of a scalar 
parameter o in order to obtain a family of controllers K(s,o) where o; variance of the 
fictitious noise injected into the input of the plant, represents plant uncertainty, and plays 
the role of a tuning parameter of the optimal observer. 


It can be shown that Loop Transfer Recovery (LTR) at the control input is achieved 


if [Ref. 6] 


K(s,o)P(s) > Li(s) (62) 


pointwise in s as 6 > o. 

The exact design procedure depends on the point where the unstructured 
uncertainties are modeled and where the loop is broken to evaluate the open-loop 
transfer matrices. Commonly either the input point or the output point of the plant is 


taken as such a point. The design which was proposed in Chapter III is shown with the 


loop broken at the control input in Figure 17. 








ACCNOISE GYRONOISE 






X= Ax + Bu + EW 
= Cx +Du + V 


Missile Model 







Optimal 
Feedback 
Gain 


Cmdacc 


Figure 17. Optimal Linear Output Feedback Control System with Loop Broken at the 


Control Input. 





4] 








It is important to realize that by adding an observer to the control system, the 


transfer function from the reference input to the reference output is not affected if the 
loop is broken at the control input. However, if the loop is broken anywhere else, the 
observer will affect the transfer function. In order to improve the robust stability of the 
LQG design proposed in Chapter III, loop transfer recovery will be applied by breaking 


the loop at the uncertainty, as shown in Figure 18. 








Optimal Kalman Filter 
Feedback 
Gain " N State Command 
xX |Matrix 
Cmdacc 


Figure 18. Optimal Linear Output Feedback Control System with Feedback Perturbation 


and Loop Broken at the uncertainty. 
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To make the design less sensitive to the perturbation A , fictitious noise will be 
added to the system at the uncertainty in the design process. Actually, fictitious noise is 
considered only in obtaining the optimal estimator gain, and is not present in the plant. 
The configuration shown in Figure 18 recalls the general A-P-K framework presented in 


Chapter II and illustrates this technique (see Figure 19). 


Fictitious Noise 


Kis) 


Figure 19. General A-P-K Framework with Fictitious Noise Added at the Uncertainty. 





In this case the assumption of fictitious white noise, in addition to any white noise 
that may actually be present, enhances the robustness as the intensity of the fictitious 


noise tends to infinity asymptotically. 








As the plant copes with the fictitious noise, the LQG controller becomes more 
robust to gain and phase plant input, but the optimality for the original nominal stochastic 
model is no longer guaranteed. In the case of non-minimum phase plants,o — «x will not 
lead to loop recovery. Thus, to recover open loop transfer functions with a series 
compensator, a plant must be minimum phase. 

The following example includes a servo in parallel with the vibrational mode 
designed with LQG techniques. The sequence illustrates how the good robustness 


properties of an LQR controller (Gain Margin = 2, Phase Margin > 60°, 6 dB margin 


against gain reductions) can be recovered by injecting fictitious noise at the control input. 








Gain Margin (GM)=1.11 
Down-Side Gain Margin (OSGM)=0.72 
Phase Margin (PM)=2.5 


Figure 20. Nyquist Polar Plot for Optimal Control. 
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Variance of Fictitious Noise = 0.001 
GM =3.5 

DSGM = .66 

PM = 34.5 


Figure 21. Nyquist Polar Plot for Optimal Control Plus Kalman Filter. 











Variance of Fictitious Noise = 0.008 
GM = 20 

DSGM = .34 

PM = 67.5 


Figure 22. Nyquist Polar Plot for Optimal Control Plus Kalman Filter. 
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Variance of Fictitious Noise = 10 
GM = 0 

DSGM = 0 

PM = 143 


Figure 23. Nyquist Polar Plot for Optimal Control Plus Kalman Filter. 





In the case of the Linear Quadratic Gaussian Controller Design for the flexible 
missile the uncertainties were not taken into account. As a result, Nominal Performance 


was achieved but Robust Performance and Robust stability were not. In other words, the 
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LQG Design becomes sensitive to parameter variation or uncertainties. The acceleration 


time response and the results of the y—analysis test are presented in Figure 24. 





0 0.5 1 
Time (sec) Freq (Hz) 
a. b. 


Figure 24. Acceleration Time Response (a) and p:-analysis (b) for LQG Design without 


Uncertainties. 





The performance of the LQG design in the presence of fictitious noise injected at 
one or two uncertainties was evaluated. Loop Transfer Recovery through injection of 
fictitious noise at M, in the presence of an uncertainty at Z, 1s applied to enhance the 
robustness of the design. The effect of this technique on the time response and resulting 
p-analysis tests are illustrated in Figures 25 through 28. 

It can be seen in Figure 25(a) that once an uncertainty at Z, is considered, the LQG 
controller is unable to track the desired reference input. In addition, robust performance, 


robust stability, and nominal performance requirements are not met (Figure 25(b)). 
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“0 0.5 1 
Time (sec) Freq (Hz) 


a. b. 


Figure 25. Acceleration Time Response (a) and p:-analysis (b) for LQG Design in the 


Presence of an Uncertainty in Z,. 


As shown in Figure 26, time response is improved through the addition of fictitious 
noise at M,, and p-analysis results indicate a corresponding improvement in robust 


stability. Nominal performance and robust performance, however, are not met. 


Amplitude 





0.5 
Time (sec) Freq (Hz) 


a. b. 


Figure 26. Acceleration Time Response (a) and j:-analysis (b) for LQGLTR Design in the 


Presence of an Uncertainty in Z, (Variance of the Fictitious Noise =1). 





50 











As the variance of the fictitious noise added at M, increases, the acceleration time 
response of the system tends to track the desired reference input, and robust stability is 


achieved, as shown in Figures 27 and 28. Once again, neither robust performance nor 


nominal performance are met. 


Amplitude 








or 0.5 1 
Time (sec) Freq (Hz) 
a. b. 


Figure 27. Acceleration Time Response (a) and :-analysis (b) for LQGLTR Design in the 


Presence of an Uncertainty in Z, (Variance of the Fictitious Noise=100). 
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Amplitude 





“0 0.5 1 


Time (sec) Freq (Hz) 
a. b. 


Figure 28. Acceleration Time Response (a) and ,.-analysis (b) for LQGLTR Design in the 


Presence of an Uncertainty in Z, (Variance of the Fictitious Noise=500). 





So far, robust stability is achieved with enough intensity of the fictitious noise but 
as a consequence of loosing nominal performance. It must be recognized that there exists 
a tradeoff between performance loss for the nominal model and robustness gain. The 
simultaneous achievement of robust performance, robust stability, and nominal 
performance will require application of more sophisticated techniques, such as linear 
quadratic gaussian loop transfer recovery with formal synthesis and H,, optimal control, 


and is beyond the scope of the present work. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

This thesis presents the design of a controller for a flexible missile model using 
Linear Quadratic Gaussian techniques in combination with Loop Transfer Recovery in 
order to enhance the Robust Stability of the design. 

A linearized model of the tail controlled surface-to-air missile was provided. The 
design of the Linear Quadratic Gaussian controller was based on a combination of the 
solution of the Stochastical Linear Quadratic Regulator (LQR) and Optimal Estimator 
problem (Kalman Filter). In order to transform the problem into a tracking problem a 
State Command Matrix ( N,) was designed to drive the reference input to a reference 
state. The Closed Loop Optimal Output Feedback Control System was then assembled 
and Robust Performance, Robust Stability, and Nominal Performance of the system was 
evaluated. 

A Linear Quadratic Gaussian controller with Loop Transfer Recovery Techniques 
was then designed to improve the Robust Stability of the design. The u—Analysis and 
Synthesis Toolbox and the Control Toolbox of MATLAB were used for the assembling, 


design, analysis and simulation of the missile system. 
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B. CONCLUSIONS 

In this study, a Linear Quadratic Gausian Controller is developed for the model of a 
tail controlled surface-to-air missile without taking in consideration the uncertainties 
present in the model. All of the required Time Domain Specifications for the missile 
autopilot are met by the LQG design. The controller does not saturate the tail deflection 
actuator rate capabilities nor destabilize the high frequency flexible body modes of the 
missile. Nominal Performance is met by the LQG design, but Robust Performance and 
Robust Stability are not met. In short, the LQG controller meets the design goal at the 
nominal operating condition but it becomes sensitive with parameter variation. Some 
linear combination of the uncertainties might destabilize the system. 

To improve the Robust Stability of the system, Loop Transfer Recovery is then 
applied to the LQG design by the injection of fictitious noise at different uncertainties. 
All of the required Time Domain Specifications for the missile autopilot are met by the 
LQGLTR design. The controller does not saturate the tail deflection seals rate 
capabilities nor destabilize the high frequency flexible body modes of the missile. Robust 
Stability is also met but as a consequence of loosing Nominal Performance. Robust 
Performance is not met. A trade-off must be observed between the loss of nominal 


Performance and the gain of Stability Robustness. 
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C. RECOMMENDATIONS 

The LQG controller designed in this work is not recommended for the missile 
model in question since it is sensitive under parameter variations or uncertainties which 
could destabilize the system. 

The LQGLTR controller designed in this work can be used for the missile model in 
question if Robust Stability is desired and the designer is willing to make a trade-off in 
the nominal performance of the system. 

The use of more sophisticated techniques such as H, Optimal control and 
u-synthesis may allow the designer to meet the design specifications as well as all of the 


performance tests. 
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